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Abstract 

The facultative anaerobe, Bacillus cereus, causes diarrheal diseases in humans. Its ability to deal with oxygen availability is 
recognized to be critical for pathogenesis. The B. cereus genome comprises a gene encoding a protein with high similarities 
to the redox regulator, Rex, which is a central regulator of anaerobic metabolism in Bacillus subtills and other Gram-positive 
bacteria. Here, we showed that B. cereus rex is monocistronic and down-regulated in the absence of oxygen. The protein 
encoded by rex is an authentic Rex transcriptional factor since its DNA binding activity depends on the NADH/NAD^ ratio. 
Rex deletion compromised the ability of 6. cereus to cope with external oxidative stress under anaerobiosis while increasing 
6. cereus resistance against such stress under aerobiosis. The deletion of rex affects anaerobic fermentative and aerobic 
respiratory metabolism of 6. cereus by decreasing and increasing, respectively, the carbon flux through the NADH-recycling 
lactate pathway. We compared both the cellular proteome and exoproteome of the wild-type and Arex cells using a high 
throughput shotgun label-free quantitation approach and identified proteins that are under control of Rex-mediated 
regulation. Proteomics data have been deposited to the ProteomeXchange with identifier PXD000886. The data suggest 
that Rex regulates both the cross-talk between metabolic pathways that produce NADH and NADPH and toxinogenesis, 
especially in oxic conditions. 
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Introduction 

Bacillus cereus is a Gram-positive, facultative-anaerobe, rod- 
shaped endospore-forming tiuman pathogen. Most of the reported 
illnesses involving B. cereus are food-borne intoxications, classified 
as emetic and diarrheal syndromes [1,2]. Diarrheal disease is due 
to vegetative outgrowth and secretion of various extracellular 
factors, including enterotoxins [3]. The most extensively studied 
diarrheal enterotoxins are hemolysin BL (Hbl), nonhemolytic 
enterotoxin (Nhe), and cytotoxin K (CytK) [4,5]. These entero- 
toxins are secreted via the Sec translocation pathway [6] . Although 
Hbl, Nhe and CytK are currently considered as the etiologic 
agents of diarrheal syndrome, other toxins, such as EntA, EntB 
and EntC, may also contribute to the pathogenicity of B. cereus 
[7,8]. To grow and produce virulence factors in the human 
intestine, B. cereus must adapt its metabolism, and regulates its 
proteome [9] in response to changes in oxygen availability. 
Indeed, B. cereus encounters oxic conditions in zones adjacent to 
the mucosal surface [10] and anoxic condition in the intestinal 
lumen [11]. Changes in oxygen availability can influence the 
relative levels of the dinucleotide, NAD'*^ and NADH, in the cell, 
and such changes are sensed by the transcriptional regulator. Rex, 
in B. suhtilis [12] as in other Gram-positive bacteria 
[13,14,15,16,17,18,19]. Depending on the cellular NADVNADH 



ratio. Rex regulators modulate the expression of genes involved in 
fermentative metabolism, biofilm formation, and oxidative stress 
[18,20,21]. Structural studies of Rex proteins have identified 
dinucleotide-binding pockets in the C-terminal domain of the 
protein. NADH binding in this region leads to a conformational 
change in the Rex homodimer, triggering a displacement of Rex 
from its recognition sites on DNA, and thus leading to de- 
repression of the downstream genes [15,19]. A Rex homologue 
has been detected in the cellular proteome of B. cereus [9]. 
Furthermore, we found a canonical Rex binding motif [15] 
overlapping ResD and Fnr binding motifs [22,23,24] in the IdhA 
promoter region [25]. In B. cereus, ResD, Fnr and LdliA regulate 
both catabohsm and enterotoxin production under both aerobic 
and anaerobic growth conditions, probably through a regulatory 
complex [26]. Rex could be, thus, a regulator of both catabohsm 
and toxinogenesis in B. cereus. A role of Rex in toxinogenesis has 
not yet been reported in B. cereus or in any other organism. 

In this study, we show that Rex from B. cereus, like its 
orthologues, is a transcriptional factor capable of interacting with 
DNA in an NADH/NAD''"-responsive manner. We demonstrate 
that B. cereus Rex is a key regulator of anaerobic fermentation, 
aerobic respiration, resistance against external reactive oxygen 
species, and toxinogenesis, by modulating the cellular and 
extracellular proteome in an oxygen-dependent manner. This 
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study provides the most comprehensive experimental information 
on proteins whose synthesis was changed in the presence of Rex. 
All together, our results offer new information about the metabolic 
events that maximize B. cereus growth in environments with 

varying oxygen conditions. 

Materials and Methods 

Bacterial strains, media, and growth conditions 

Escherichia coli TOP 10 (Invitrogen) was used as the host for 
cloning experiments, and E. coli SCSllO (Stratagene, La JoUa, 
CA) was used to prepare DNA for B. cereus transformation. B. 
cereus ATCC 14579 [27] was used as the parent strain for the 
construction of the rex deletion mutant. E. coli strains were grown 
at 37°C, with agitation, in Luria broth (LB). B. cereus strains were 
cultured in batches (three independent cultivations per strain) at 
two oxygen availabilities, i.e. pO2 = 0% and pOj = 100% [8,9]. 
Each batch culture was inoculated with a subculturi; grown 
overnight at an initial optical density at 560 nm (ODsgo) equal to 
0.02. For B. cereus cultivation, the minimal MOD medium was 
supplemented with 30 mM glucose as the carbon source [28]. 
Anaerobic and aerobic batch cultures were performed at 37°C in a 
2 L bioreactor (BioFlo/CcUiGen 115, New Brunswick), and the 
working volume was maintained at 1 .8 L. The pH was kept at a 
controlled value of 7.2 by automatic addition of 5 M KOH. B. 
cereus growth was monitored spectrophotometrically at 560 nm 
and calibrated with cell dry-weight measurements as previously 
described [25]. B. cereus cells were harvested by centrifugation 
when they reached their maximal growth rate {[i = |J.,-nax) and 
immediately frozen until proteomic analysis. Supcrnatants were 
kept at — 20°C for glucose and glucose-by-product assays and 
exoproteomic analysis. 

Construction of the B. cereus Arex mutant strain 

The deletion mutant, ATCC 14579 Arex, was constructed as 
follows. Two DNA fragments encompassing the 5'untranslated 
region (UTR) and 3'UTR of rex (BC 0291) were generated by 
PGR using primer pairs, rexlF (5'-GCCATGTTAATGTTTC- 
GATGTGT-3') and rex2R (5'- GGCGGG ATGTTTTAG- 
GAGTGGCTTGTGG-3', Smal restriction site is underlined), 
and rex3F (5'- CCCGGGG TTTACTTTTTGAAAAACTATC- 
CACAA-3', Smal restriction site is underlined) and rex4R (5'- 
TGCATTATGTATCGTGCTTTGG-3'), respectively. The re- 
sulting 840 bp ySmal and 827 bp 5'SmaI DNA fragments were 
cloned into the TA cloning vector, pCR4-TOPO (Invitrogen, La 
JoUa, CA), generating plasmids, pCR4mutr(;xl and pCR4mutrex2, 
respectively. The 840 bp DNA fragment encompassing the 
5 'UTR region of rex was isolated from pCR4mutrexl with Pstl 
and Smal, and subcloned into pCR4mutr'«x2 to generate 
pCR4mutrex3. A 1.5 kbp Smal fragment containing the entire 
spectinomycin gene, spc [29], was purified from pDIA [25]. This 
purified fragment was ligated into Smal-digested pCR4mutrex3. 
The resulting plasmid, pCR4mutre«;4, was digested with EcoRl 
and the resulting 3067 bp 5'rexUTR-5^c-3'rexUTR was subse- 
quently inserted into the EcoRl site of pMAD [30]. The resulting 
plasmid was introduced into B. cereus cells by electroporation. 
The rex ORF was deleted and replaced with spc via a double- 
crossover event [30]. Chromosomal allele exchanges were 
confirmed by PCR with oligonucleotide primers located upstream 
and downstream of the DNA regions used for allelic exchange. To 
complement the rex gene in trans, a DNA fragment of 932 bp 
encompassing the rex ORF (630 bp) and its promoter region 
(219 bp) was first PCR amplified using the primer pairs, rex- 
compF (5 '- GGATCC CGTTCGAAAGCGCGTTTACTTG-3 '; 



the BamHl restriction site is underlined) and rexcompR (5'- 
GAGCTC GATTTTAATTTGGCACTTCGCC-3'; die Sad re- 
striction site is underhned), and then cloned into the pCRXL- 
TOPO plasmid (Invitrogen). The PCR fragment was then cut with 
BamHl and Sad and ligated into pHT304 [31], digested with the 
same restriction enzymes. The integrity of the recombinant vector 
(pHT304r«x) insert was verified by sequencing. 

Phenotypic characterization of B. cereus Arex using the 
API-50CHB testsystem 

The carbohydrate metabolism of the wild-type strain (WT) 
transformed or not with pHT304, Arex mutant transformed or not 
with pHT304 and Arex mutant transformed with pHT304rex was 
examined using API 50 CHB strips (BioMerieux SA, France). The 
results showed that WT and WT(pHT304) did not ferment 
turanose. In contrast, Arex and Arex (pHT304) showed a typical 
positive reaction in turanose test fermentation. Transformation of 
Arex mutant with pHT304rex inhibited the capacity of Arex to 
ferment turanose. 

Measurement of glucose and by-product concentrations 

Enzymatic test kits from Diffchamb (Lyon, France), R- 
Biopharm (Saint-Didier au Mont-d'Or, France), and Roche 

(Meylan, France) were used to analyze the glucose, lactate, 
ethanol, formate, acetate, and succinate concentrations in the 
supcrnatants of 4 mL cell cultures obtained after centrifugation at 
1 0,000 xg for 5 min (4°C). The specific glucose consumption rate, 
defined as the differential change in glucose concentration with 
time, was calculated from the equation, 5giu<:os<: = V-^Yx, where |J, is 
the specific growth rate (h ') and is the biomass yield (g.mol 
carbon substrate '). 

Gene expression analysis by RT-PCR 

RT-PCR was performed using SYBR Green technolog\- on a 
Lightcycler instrument (Roche applied Science) as described 
j)rc\i()usly [32]. The primers used in this study have been 
described previously [8]. 

Purification of Rex 

The rex ORF was amplified by PCR from B. cereus ATCC 
14579 using the oligonucleotides, pETlOlre.xF (5'-ACCATG- 
GATCAGCAAAAGATTCCA-3') and pETlOlrexR (5'- 

was introduced as a blunt-end PCR product into pETlOl/D- 
TOPO (Invitrogen). The integrity of the inserted sequence was 
confirmed by DNA scqiu-ncing. The resulting construct was 
transformed into E. coli BL21-CodonPlus(DE3)-RIL strain 
(Stratagene) for protein production. BL21-CodonPlus(DE3)-RIL 
cells carrying the pETlOl-rex expression plasmid were grown in 
1 L LB medium containing 100 |J.g.mL ' ampicUlin, at 30°C with 
agitation (200 rpm) until the cell density reached an OD^oo of 
about 0.5. After this, 0.5 mM isopropyl-P-D-thiogalactopyrano- 
side (IPTG) was added to the culture, and growth was continued 
for four hours. The cells were harvested by centrifugation 
(6500 rpm for 10 min at 4°C), washed, and resuspended in 
20 mL ice-cold extraction buffer consisting of 50 mM TRIS 
buffered at pH 8.0, 0.3 M NaCl, 10% glycerol, and a protease 
inhibitor cocktail (one tablet of Complete Mini, EDTA free, 
Roche). Cells were then incubated with 0.5 mg/mL lysozyme for 
45 min under gentie agitation at 4°C and disrupted by sonication 
for 5 min at 4°C using a Vibra Cell ultrasonicator (Fisher Bioblock 
Scientific). The lysates were clarified by centrifugation at 
9000 rpm for 30 min at 4°C and loaded onto a 2 mL 
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Co """-immobilized TALON metal affinity chromatography col- 
umn (Clontech) equilibrated with the extraction buffer. The 
column was washed with 30 mL of extraction buflFer containing 
20 mM imidazole and Rex was eluted with 3 mL extraction buffer 
containing 200 mM imidazole. The eluted fraction was desalted 
by dialysis, concentrated using Nanosep 30 kDa-molecular-mass 
cutoff devices (Omega disc membrane; Pall filtron), and stored at 
— 80°C until analysis. The purity of the protein was estimated to 
be above 90% by Coomassie blue-stained SDS-PAGE. The 
protein concentration was determined using a Bradford assay 
(Interchim) with bovine serum albumin as the reference. 

Electrophoretic mobility shift assays (EMSAs) 

Nucleic acid fragments containing the promoter regions otldhA 
and rex were PCR-amplified using biotinylated forward primers, 
LdhAF (5'-ACCTGCTAATCCGATGATTG-3') and RexF (5'- 
CAAGAATCGTTTCTGCACCG-3'), and nonbiotinylated re- 
verse primers, LdhAR (5'-GGATCCAACTAATCCAGTAC-3') 
and RexR (5'-GCTTACCAGAAAGAGATAAG-3'). The DNA 
used as negative control was a fragment of the ssuRNA BC0007 
sequence (NC_004722), which was amplified with the biotinylated 
ssubioF (5'-GGTAGTCCACGCCGTAAACG-3') and ssuR (5- 
GACAACCATGCACCACCTG-3') primer pair. The ,5'-labeled 
ampKcons were purified using the High Pure PGR Product 
Purification Kit (Roche). Binding reactions were performed for 
30 min at 37°C by incubating biotin-labeled DNA fragments 
(2 nM per reaction) with difierent amounts of Rex in 10 mM Tris- 
HCl buffered at pH 7.5, and containing 50 mM KCl, 2.5% 
glycerol, 5 mM MgCl2 and 5 mg/L poly(dI — dC). The samples 
were resolved by electrophoresis on a 6"/o nondenaturing 
polyacrylamide gel and electrotransferred onto Hybond N+ Nylon 
membranes (Amersham). Biotin-labeled DNA was detected using 
the LightShift Chemiluminescent EMSA Kit (Pierce). 

Proteomic sample preparation and nanoLC-MS/MS 
analysis of tryptic peptides 

Three independent biological replicates were harvested for each 
of the two conditions (aerobiosis and anaerobiosis) and two strains 
(Arex mutant and its parent strain, ATCC 14579). The 
extracellular proteins of the 12 cultures were extracted by 
trichloroacetic acid precipitation [8]. The cellular proteins from 
the 1 2 samples were obtained as previously described [9] . The 24 
resulting samples were subjected to SDS-PAGE, and then 
identified after trypsin proteolysis by nanoLC-MS/MS tandem 
mass spectrometry with an LTQ;Orbitrap XL mass spectrometer 
as previously described [8,9]. A total of 131 nanoLC-MS/MS runs 
were carried out to acquire the whole dataset. The MS/MS 
spectra were assigned to tryptic peptide sequences with the Mascot 
Daemon software (version 2.3.2; Matrix Science) with mass 
tolerances of 5 ppm on the parent ion and 0.5 Da on the MS/ 
MS, fixed modification for carbamidomethylated cysteine, and 
variable modification for methionine oxidation. Mascot results 
were parsed with a ^-value threshold below 0.05 for peptide 
identification and proteins were validated when at least two 
peptides were detected. The number of MS/MS spectra per 
protein recorded by nanoLC-MS/MS was extracted for each 
sample. In each condition, proteins were further considered for 
comparison only if peptides were seen in at least two of the three 
replicates. The resulting datasets were normalized taking into 
account the total protein concentration of the corresponding pellet 
and supernatant. Protein concentrations in B. cereus were 
determined using the Reducing agent Compatible Detergent 
Compatible (RCDC) protein assay (Bio-Rad) following the 
supplier's instructions. MS/MS spectral counts were compared 



with the TFold method using the PatternLab software program 
2.0.0.13 [33] using a^-value cut-offset at 0.05. Log2 (fold-change) 
were then calculated for comparisons and only the proteins for 
which thep-value was lower than the 0.05 cut-off were considered. 
The mass spectrometry proteomics data have been de- 
posited to the Pr<)t(X)m(;Xchange Consortium (http:/ /www. 
proteomexchange.org) via the PRIDE partner repository [34] 
with the dataset identifier PXD000856 and DOI 10.6019/ 
PXD000856. 

Exposure of bacteria to H2O2 and viability assays 

Oxidative stress resistance of the B. cereus Srex mutant and the 
parent strain, ATCC 14579, was assessed by exposing aerobically 
grown (OD~0.4) and anaerobicaUy grown (OD~0.2) cells to 20 
and 5 mM H2O2, respectively. Samples were taken prior to 
oxidative stress (time zero) and after 20 min. Ahquots (100 |J.l) of 
the samples were diluted in H2O, appropriate dilutions of the 
culture were plated onto LB agar, and after overnight incubation 
at 37°C the colony forming units (CFUs) were counted. AH the 
experiments were performed at least in triplicate, and at least two 
technical rephcates from each dilution step were carried out to 
determine the number of CFUs. 

Results 

Expression analysis of rex in 6. cereus ATCC 14579 at two 
different PO2 

We identified BC_0291 in the B. cereus ATCC 14579 genome 
[27] as the homologue to the Bacillus subtilis transcriptional 
repressor, Rex [12,35,36], with 90% sequence identities. We 
established by 5'RACE PGR (Fig. SI) that a transcriptional start 
site (G) was located 41 bp from the translational start site (ATG) of 
the BC_0291 open reading frame. Upstream of this transcriptional 
start site, we identified a potential housekeeping aA-type 
promoter, TATAC.-\N( 1 7)TAAAGT. The stop codon (TAA) 
overlapped an inverted repeat (AAAACGCAGAGG(N6)C- 
CTCTGCGTTTT; AG = -23.2 kcal/mol) that may be a tran- 
scriptional terminator, suggesting that rex was monocistronicaUy 
transcribed. Using real-time RT-PCR, we investigated the 
expression of B. cereus rex uiitk-r aerobiosis and anaerobiosis in 
early-growing ATCC 14579 cells. The results indicated that rex 
expression was significantly higher under aerobiosis than under 
anaerobiosis (log2 fold-change = 1.1, ^-value<0.05). 

Rex regulates glucose catabolism in 6. cereus 

Given that Rex orthologues function as NAD^/NADH- 
responsive regulators of metabolism in several Gram-positive 
bacteria (such as B. suhlilis, S. aureus and S. coelicolor 
[12,13,14,15,36,37]), we sought to determine if B. cereus rex 
contributed to anaerobic fermentative and aerobic respiratory 
growth. As shown in Table 1, rex deletion slightly increased the B. 
cereus growth rate under both anaerobiosis and aerobiosis without 
a significant change in final biomass. Under anaerobiosis, the 
increase of growth rate was not related to glucose consumption, 
suggesting perturbation of the fermentative pathways. In accor- 
dance with this hypothesis, the spectra of fermentation end 
products were significantly modified by the rex deletion (Table 1). 
We observed that formate, acetate and ethanol production were 
promoted at the expense of lactate production. In addition, the 
ethanol-to-acetate ratio increased in Arex mutant, while NADH 
recovery levels remained unchanged and ATP yields tended to 
increase. Under aerobiosis, respirator)' Arex cells secreted ~8-fold 
higher amounts of lactate compared with the wild-type cells, while 
producing lower levels of acetate. We conclude that rex deletion 
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affects the carbon flux through the NADH recycUng lactate 
pathway under both anaerobiosis and aerobiosis (Fig. 1). Com- 
plementation for gTowth of the non-polar rex mutant in B. cereus 
ATCC 14579 was not possible using a pHT304-based plasmid. 
This is probably due to the presence of multiple copies of the rex 
promoter region and/or overexpression of rex [38]. Therefore, to 
validate the role of Rex in aerobic respiration and anaerobic 
fermentation in B. cereus, we deleted the rex gene of B cereus 
F4430/73 [39] and cultured the mutant strain under both 
aerobiosis and anaerobiosis. The results showed that the growth 
phenotype of Arex was the same in strains ATCC 14579 and 
F4430/73 (Table SI). 

Rex regulates hydrogen peroxide resistance of B. cereus 
cells 

To evaluate the impact of Rex deficiency on B. cereus resistance 
to oxidative stress, aerobicaUy and anaerobicaUy grown ATCC 
14579 cells were exposed to hydrogen peroxide (20 and 5 mM 
H2O2, respectively). Figure 2 shows that aerobicaUy grown Arex 
cells were less susceptible to H2O2 harmful effects than WT cells, 
while the anaerobicaUy grown Arex ceUs were more susceptible. 
These data indicate that Rex restricts the resistance to oxidative 
stress of aerobicaUy grown B. cereus cells, whUe sustaining a high 
resistance to oxidative stress of anaerobicaUy growing cells. To 
validate the role of Rex in the B. cereus resistance to external 
H202_ we repeated the experiments using the F4430/73 Arex 
strain (Fig. S2). The results showed that Rex deficiency simUarly 
impacts F4430/73 and ATCC 14579 ceUs whUe the survival of 
F4430/73 anaerobicaUy grown cells was strongly higher than the 
survival of ATCC 14579 ceUs. 



DNA binding activity of 6. cereus Rex 

Previous studies have shown that proteins belonging to the Rex 
famUy function as dimers that bind to promoter regions containing 
DNA motifs with the 5'-TTGTGAAnnnnTTCACAA-3' consen- 
sus sequence [21]. This typical binding motif is missing in the rex 
regulatory region of 5. cereus, as reported for S. aureus and B. 
suhtilis [12,15]. However, a putative binding motif with two 
mismatches compared with the known consensus motif was found 
upstream of the transcription start site ofldhA [25] in 5. cereus, as 
in S. aureus and B. suhtilis. To test whether B. cereus Rex binds to 
the IdhA promoter region, we overexpressed B. cereus rex in E. 
coli, purified the dimeric His6-tagged recombinant protein (Fig. 
S3), and performed electrophoretic mobUity shift assays (EMSA). 
As a negative control, we assessed the binding of the Rex protein 
to the ssu DNA fragment (Fig. 3A). The results showed that a 
concentration of 0.6 |iM Rex led to a complete shift of the IdhA 
fragment (2 nM). As expected, no shift was observed with the rex 
DNA fragment. The specificity of the binding was demonstrated 
by the absence of DNA shift with the negative control. The second 
EMSA experiment (Fig. 3B) showed that 10 mM NADH impaired 
Rex binding to the IdhA promoter region while 100 mM NAD"*" 
did not interfere with Rex-DNA complex formation. An NAD"""/ 
NADH ratio of five led to an incomplete shift of the IdhA 
fragment. We conclude that in B. cereus, as in B. suhtilis and S. 
aureus [15], Rex binding activity depends on the NADH/NAD"^ 
ratio. 

Global comparative proteomics of 6. cereus Arex and 
wild-type cells 

To investigate the regulatory role of Rex in aerobic respiratory 
and anaerobic fermentative growth of B. cereus and decipher its 
role in toxinogenesis, we carried out a comparative proteome 



Table 1. Results from controlled batch cultures of Arex mutant and Its parent strain, 6. cereus ATCC 14579^. 







Anaerobic fermentative growth 


Aerobic respiratory growth 




wt" 


Arex 


wt" 




Maximal specific growth rate (fimax) (h^^) 


0.84 ±0.04 


0.90±0.02» 


1.22±0.11 


1.41 ±0.04* 


Final biomass (g. liter"') 


0.78±0.01 


0.83±0.06 


2.38±0.12 


2.31 ±0.12 


Ygiucose (g of cells, mol of glucose"') 


26±1 


28±1 


79±1 


77±1 


Maximal specific glucose consumption(mmol.g"Vh"') 


32±2 


32±2 


15±1 


18±2 


Yields of end products (mol. mol glucose 


Lactate (Vi^MtJ 


1.33±0.01 


1.25±0.02* 


0.05±0.02 


0.40±0.01* 


Acetate (Y acetate) 


0.32±0.01 


0.36±0.01* 


1.05±0.01 


0.69±0.03* 


Formate (Yformate) 


0.23±0.01 


0.27±0.01* 


0.01 ±0.01 


0.01 ±0.01 


Ethanol (Yethanoll 


0.13±0.01 


0.18±0.01* 


NZ" 


NZ 


Succinate (/;„„;„„„) 


0.01 ±0.01 


0.01 ±0.01 


0.05 ±0.01 


0.04±0.01 


Ethanol versus Acetate 


0.40 


0.50* 






ATP yield" 


2.10 


2.15 


nd' 


ND 


NADH recovered^ 


1.2 


1.2 


ND 


ND 



^Cells were grown under anaerobiosis (p02 = 0%) and full aerobiosis (p02 = 1 00%). Data are the means of triplicate measures obtained from three independent cultures. 

^WT, wild-type parent strain B. cereus ATCC 14579. 

'^Yields of end products were calculated at the stationary phase. 

■^NZ, yield was below 0.01 mol. mol glucose~\ 

^ATP yield was calculated as moles of ATP produced per mole of consumed glucose, and was equal to Y\actate+yethano&^*y acetate- 
'^ND, not determined. 

^NADH recovery was calculated as the ratio of pathways producing NADH versus those consuming NADH (producing NAD^), and was equal to {lactate+2xacetate+ 
2 X ethanol —formate)/(lactate+2 xethanol). 
*p<0.05 versus WT in Student's t-test. 
doi:l 0.1 371/journal.pone.01 07354.t001 
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PTS transport system 



D-Glucose 



Phosphoenol 
pyruvate • 



Glycolysis 



Butane-2,3-diol9-H; 




PDHc 



BC4S1G I 
BC4517 I 

Ac*^ ©Succinate 
^ — BCSSS— I 



Figure 1. An overview of the main anaerobic and aerobic glucose catabolic pathways utilized by B. cereus ATCC 14579. The proteins 
detected in tinis study are indicated by their BC number. Protein names and functions are listed in Table S6. The form filling indicates fold-change 
values that satisfied the Student's f-test statistical criteria (p-value<0.05) in anaerobiosis (blue) and aerobiosis (yellow). Red and green BC numbers 
indicated significant increase and decrease, respectively, of abundance level of the proteins in Arex mutant compared with wild-type. The 
nicotinamide nucleotides are indicated in red (NADP^/NADPH) or blue (NADVNADH). 
doi:10.1371/journal.pone.0107354.g001 



analysis from cells harvested in exponential growth phase 
(|I = |J,„iax)- We compared the cellular and extracellular proteomes 
of aerobically and anaerobicaUy grown Arex with those of its 
parental strain, with biological triplicates for each condition. A 
total of 1,620 proteins were identified in the whole study (Tables 
S2, S3, S4 and S5), with 1,450 proteins specifically detected in the 
whole-cell shotgun analysis and 170 further proteins observed in 
the exoproteome study. These proteins were quantified by spectral 
count as previously described [8,9]. The Rex-dependent changes 
in protein abundance, expressed as loga of the ratio of a protein's 
abundance at a given condition (aerobiosis or anaerobiosis) 



relative to the wild-type, are presented in Tables S6 and S7, 
where the proteins with significant abundance changes based on 
the PatternLab Tfold student t-test {p-value below 0.05) are 
highlighted while proteins with non-significant changes are in grey 
color. 

(i) Anaerobiosis. Among the whole cellular proteome (Table 
S6), 145 proteins showed significant abundance level changes in 
Arex cells relative to WT (^-values<0.05): 89 proteins were 
significandy up-regulated and 56 proteins were significantly down- 
regulated in Arex cells. In the exoproteome (Table S7), 41 proteins 
showed significant abundance level changes in Arex cells relative 
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Figure 2. Survival of B. cereus ceWs towards external hydrogen peroxide under aerobiosis and anaerobiosis. Cells were grown In liquid 
cultures to the mid-exponentlal growth phase either aeroblcally or anaeroblcally, and subjected to 20 or 5 mM H2O2, respectively. Samples were 
taken at 0 and 20 min after H2O2 addition. Colony forming units per mL were counted and expressed as (N/No)x100. Error bars represented the 
standard deviation from three Independent measures. Significant differences (p-value<0.05) between mutant and wild-type strains are Indicated 
with asterisks. 

dol:1 0.1 371 /journal.pone.01 07354.g002 



to WT: 1 7 proteins were significantly up-regulated and 24 proteins 
were significantly down-regulated in Arex cells. 

(ii) Aerobiosis. Table S6 indicates that 132 proteins showed 
significant abundance level changes in Arex cellular proteome (p- 
values<0.05): 82 proteins were significantly up-regulated and 50 
proteins were significantly down-regulated in Arex cells. Table S7 
indicates that 70 proteins showed significant abundance level 
changes in Arex cellular exoproteome (|!)-values<0.05): 35 were 
significantly up-regulated and 35 were significantly down-regulat- 
ed in Arex cells. 



(iii) Global analysis (Fig S4). Only 8 cellular proteins 
showed similar behavior under both aerobiosis and anaerobiosis (4 
were up-regulated and 4 were down-regialated in Arex cells). Only 
3 extracellular proteins were down-regulated under both anaero- 
biosis and aerobiosis. However, all the extracellular proteins that 
were up-regulated under anaerobiosis were also up-regulated 
under aerobiosis. Taken together, the data indicate that Rex (i) 
differently modulates the cellular proteome of B. cereus anaero- 
blcally and aeroblcally grown cells in terms of the identity of 
proteins that show significant abundance level changes, (ii) acts 
mainly as a repressor at the cellular proteome level under both 
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Figure 3. /n vitro binding of Rex to IdhA promoter region determined by EMSA. Panel A. Blotln-labeled PCR products (2 nM) 
corresponding to IdhA and rex promoter regions and a negative control (a fragment of ssuRNA BC_0007 sequence) were mixed with Increasing 
concentrations of purified Rex. Panel B. PCR product corresponding to the IdhA promoter region was Incubated with 0.6 |iiVl purified Rex and 
different concentrations of NAD* and NADH, as Indicated. 
dol:1 0.1 371/journal.pone.01 07354.g003 
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Table 2. Cellular abundance level changes of cell membrane, transport, stress response, transcriptional regulation, virulence and 
phage-related proteins in Arex mutant relative to wild-type under anaerobiosis and aerobiosis. 



Abundance level change and significance 
Anaerobiosis Aerobiosis 
Accession no Protein Gene log2 log2 



(NP) 


name 


number 


Functionnal annotation 


(fold-change)^ 


p-value 


(fold-change) 


p-value 


Cell wall, membrane and cell surface proteins 


NP_833630.1 




BC_3910 


N-acetylglucosaminyl transferase 


1,561 


0,018 


-1,515 


0,197 


NP_835080.1 


TagA 


BC_5419 


N-acetyl-beta-D-mannosaminyltransferase 


1,263 


0,042 


-0,377 


0,382 


NP_832677.1 




BC_2929 


Peptidoglycan N-acetylglucosamine deacetylase 


0,832 


0,019 


0,799 


0,183 


NP_835088.1 




BC_5427 


Cell wall biosynthesis glycosyltransferase 


1,120 


0,046 


-0,889 


0,070 


NP_834932.1 


AmsK 


BC_5269 


Amylovoran biosynthesis AmsK 


-0,713 


0,112 


1,595 


0,001 


NP_834404.1 




BC_4698 


Choline kinase 


0,949 


0,032 


0,214 


0,199 


NP_834827.1 


Lgt 


BC_5163 


Prolipoprotein diacylglyceryl transferase 


1,007 


0,012 


0,81 1 


0,004 


NP_834936.1 




BC_5273 


UDP-bacillosamine synthetase 


1,322 


0,033 


0,632 


0,110 


NP_835093.1 




BC_5432 


Bactoprenol glucosyl transferase 


1,218 


0,032 


0,536 


0,306 


NP_835019.1 




BC_5358 


Collagen adhesion protein 


0,918 


0,088 


1,170 


0,045 


NP_830682.1 




BC_0896 


S-layer protein 


-0,234 


0,180 


0,286 


0,017 


Transport-related proteins 


NP_829996.1 


SecE 


BC_0115 


Preprotein translocase subunit 


0,971 


0,044 


0,411 


0,252 


NP_834447.1 




BC_4743 


ABC transporter ATP-binding protein 


1,561 


0,018 


-0,029 


0,485 


NP_834539.1 




BC_4839 


ABC transporter ATP-binding protein 


1,585 


0,031 


0,856 


0,115 


NP_830141.1 


Uup 


BC_0290 


ABC transporter ATP-binding protein 


0,151 


0,382 


0,705 


0,045 


NP_830187.1 




BC_0348 


Methionine ABC transporter 


1,690 


0,004 


2,160 


0,064 


NP_830967.1 


OppD 


BC_1 1 82 


Oligopeptide transport ATP-binding protein 


0,651 


0,002 


-0,044 


0,421 


NP_830660.1 




BC_0874 


Arginlne ABC transporter ATP-blndIng protein 


-0,396 


0,349 


1,218 


0,024 


NP_832817.1 


Cute 


BC_3071 


Copper tansport 


0,903 


0,007 


-0,358 


0,341 


NP_830432.1 




BC_0615 


Di-/trlpeptlde transporter 


0,837 


0,047 


0,422 


0,232 


NP_830429.1 




BC_0612 


L-lactate permease 


0,748 


0,026 


-0,644 


0,261 


Stress response-related proteins 


NP_831381.1 




BC_1 603 


Cold shock protein 


1,120 


0,037 


-0,474 


0,277 


NP_830286.1 


TelA 


BC_0447 


Tellurite resistance protein 


0,084 


0,432 


0,840 


0,008 


NP_833675.1 




BC_3956 


GTP-binding protein TypA/BlpA 


0,310 


0,049 


0,043 


0,419 


NP_833972.1 




BC_4258 


Hydroxyacylglutathlone hydrolase 


0,824 


0,028 


-0,578 


0,213 


NP_833576.1 




BC_3855 


Alkaline-shock protein 


0,014 


0,488 


0,740 


0,026 


NP_834187.1 


OhrA 


BC_4475 


Organic hydroxyperoxide protein A 


0,189 


0,340 


2,029 


0,015 


NP_829959.1 


HsplS 


BC_0062 


Heat shock protein 15 


0,799 


0,082 


1,667 


0,007 


NP_830216.1 


AhpC 


BC_0377 


AlkyI hydroperoxide reductase C22 


-0,136 


0,426 


1,029 


0,007 


NP_830941.1 


KatE 


BC_1 1 55 


Catalase 


1,111 


0,014 


-0,074 


0,466 


NP_834345.1 




BC_4639 


Thiol peroxidase 


-0,044 


0,463 


0,737 


0,034 


Transcriptional Regulation 


NP_833888.1 


ArgR2 


BC_4174 


Arginlne biosynthesis repressor 


1,511 


0,030 


1,269 


0,079 


NP_833771.1 




BC_4053 


Transcriptional regulator, GntR family 


0,895 


0,019 


0,227 


0,399 


NP_833799.1 




BC_4081 


MarR family transcriptional regulator 


0,029 


0,464 


1,903 


0,001 


NP_834861.1 




BC_5197 


MarR family transcriptional regulator 


0,345 


0,029 


0,604 


0,123 


NP_832113.1 




BC_2351 


MerR family transcriptional regulator 


0,485 


0,276 


1,170 


0,045 


NP_830532.1 


TenA 


BC_0742 


Transcriptional activator 


0,791 


0,046 


-0,621 


0,193 


NP_834137.1 




BC_4425 


Transcriptional regulator 


0,632 


0,044 


-1,029 


0,051 


NP_831310.1 




BC_1531 


Transcriptional regulatory protein 


-0,029 


0,471 


1,480 


0,004 


NP_831 256.1 


ResD 


BC_1477 


Transcriptional regulatory protein 


0,356 


0,008 


-0,599 


0,119 


NP_830596.1 




BC_0806 


BigG family transcription antiterminator 


0,485 


0,276 


1,922 


0,000 
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Table 2. Cont. 





Abundance level change and significance 










Anaerobiosis 




Aerobiosis 




Accession no 
(NP) 


Protein 
name 


Gene 
number 


Functionnal annotation 


log2 

(fold-change)° 


/>value 


loga 

(fold-change) 


p-value 


Virulence-related proteins 


NP_834769.1 


Hlyl 


BC_5101 


Hemolysin 1 


-0,152 


0,438 


2,611 


0,001 


NP_832488.1 


Npr2 


BC_2735 


Bacillolysin 


-1,089 


0,193 


2,120 


0,013 


NP_833332.1 


HhoA 


BC_3600 


Protease 


0,485 


0,098 


1,690 


0,026 


NP_83 1147.1 




BC_1 366 


SSEB protein 


1,163 


0,036 


1,269 


0,055 


Phage-related proteins 


NP_831667.1 




BC_1 894 


Phage protein 


3,023 


0,014 


-1,358 


0,286 


NP_831632.1 




BC_1859 


Phage protein 


2,526 


0,039 


-0,644 


0,261 


NP_831635.1 




BC_1862 


Phage protein 


1,714 


0,027 


-0,761 


0,302 


NP_831673.1 




BC_1901 


Phage protein 


1,727 


0,012 


-1,889 


0,205 


NP_831676.1 




BC_1 904 


Phage protein 


1,975 


0,038 


-0,515 


0,367 



Only changes satisfying statistical criteria (p-value<0.05) at least in one growth condition are shown. 

^The relative amount of each protein was determined using PatternLab software. Numbers in bold indicate that satisfied statistical criteria. See Table S6 for details on all 
the proteins showing significant abundance level changes in Arex cells. 
doi:l 0.1 371/journal.pone.Ol 07354.t002 



aerobiosis and anaerobiosis and, has a stronger impact on the 
exoproteome of aerobically grown than anaerobically grown cells. 

Insights into the cellular proteome of the Arex mutant 

The list of cellular proteins showing significant abundance 
changes (j5)-values<0.05) due to the absence of Rex includes 
proteins related to glucose metabolism, amino acid, nucleotide and 
lipid metabolism, protein folding, response to stress, virulence, cell 
wall/ membrane assembly, transport, and a number of proteins of 
unknown function under both aerobiosis and anaerobiosis (see 
Table S6 for details on protein characteristics, the fold-changes 
and j&-values in each comparison). 

(i) Impact of Rex on glucose metabolism-related proteins 
under anaerobic fermentative conditions. Figure 1 shows 
an overview of the glucose catabolic pathways of B. cereui. As 
illustrated in Fig. 1, the abundance level of key fermentative 
enzymes, such L-lactate dehydrogenases (Ldh; BC_1924, 
BC_4870, BC_4996), pyruvate formate lyase (PflB; BC_0491), 
alcohol dehydrogenase (BC_4365, BC_2220, BC_0802) and 
acetate kinase (BC_4634) were not significandy different in the 
mutant compared with wild-type cells. However, the lack of Rex 
resulted in a decrease of the abundance of the IIA component of 
the glucose-specific phosphotransferase (PTS) system (BC_5320, 
log2 = —0.4; p-vahit = 0.02 1). A significant decrease in abundance 
was also observed for the pyruvate-formate-lyase-activating 
enzyme, PflA (BC_0492; log2 = -1; j!)-value = 0.048). PflA is the 
sole enzyme able to activate the oxygen-sensitive PflB (BC_0491), 
which converts pyruvate to acetyl-CoA and formate (Fig. 1) [40]. 
In contrast, anaerobically grown Arex cells sustain a higher level of 
phosphoglucomutase (Pgm; BC_4919) than the wild-type cells 
(log2 = 0.6; ^-value = 0.036) and thus, could sustain a higher 
carbon flux through gluconeogenesis than glycolysis (Fig. 1). As for 
Pgm, the malate dehydrogenase (Mdh; BC_4592) abundance level 
was increased in the absence of Rex under anaerobiosis (log2 =1.7; 
j!) -value = 0.017). Mdh is one of the enzymes involved in the 
NADH-dependent reduction of oxaloacetate to malate, which is 
the first reaction in the succinate pathway (Fig. 1). The increase of 



Mdh abundance in the Arex mutant compared with wild-type was 
associated with the increase of NADPH-producing malic enzyme 
(ME; BC_4604; log2 = 0.5; ^-value = 0.033), which converts 
malate into pyruvate. Taken together, our data suggest that Rex 
may modulate the abundance level of key metabolic enzymes to 
control carbon flow through fermentative pathways, possibly at the 
expense of the pathway that couples NADH-recycling Mdh with 
NADPH-producing ME. 

(ii) Impact of Rex on glucose metabolism-related proteins 
under aerobic respiratory conditions. Figure 1 shows that 
ME was the sole enzyme showing an abundance level change in 
both aerobically and anaerobically grown mutant cells (log2 — 0.4; 
^-value = 0.01 under aerobiosis). 6-Phosphogluconolactonase (Pgl; 
BC_3368) showed a higher increase than ME (log2=1.3; p- 
value = 0.047) in Arex aerobically grown cells. The increase of the 
Pgl abundance level may have an impact on the activity of the 
NADPH-generating pentose phosphate pathway (PPP) in the Arex 
mutant (Fig. 1). Two glycolytic proteins showed significant 
abundance decreases in aerobically grown mutant cells: the 
NADPH-producing glyceraldehyde-3-phosphate-dehydrogenase 
(NADPH-GAPDH; BC_0868; (log2= -0.5; jf) -value = 0.012) and 
pyruvate kinase (PK; BC_4599; (log2=-0.4; /;-value = 0.032), 
which catalyzes the last steps of glycolysis. The decrease of the 
NADPH-GAPDH level could increase the carbon flow through 
the NADH-producing GAPDH (BC_4583, BC_5333) and thus 
stimulate flux through the upstream glycolysis pathways, especially 
through the NADH-recycling lactate pathway (Fig. 1). Finally, 
Rex may modulate glycolytic flux by controlling the abvmdance 
level of key enzymes of glycolysis and pentose phosphate pathway 
when cells are grown under aerobiosis in glucose-containing 
medium. 

(iii) Impact of Rex on other relevant cellular 
pathways. Cell-wall and cell surface relaled-proteim. Among 
the 20 proteins classified in the cell-waU/ surface functional group, 
8 showed significant abundance level changes (p<0.05) in 
anaerobically grown Arex mutant cells while 4 proteins showed 
significant abundance level changes in aerobically grown cells 
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(Table 2). A protein annotated as a putative prolipoprotein 
diacylglyceryl transferase (Lgt; BC_5163), which may be required 
for lipid modification of the cysteine residue present within the 
lipobox of prolipoproteins [41], showed a significant change of 
abundance level under both conditions: it was more highly 
detected in anaerobically grown cells and less detected in 
aerobicaUy grown cells when Rex was absent. 

Transport-related proteins. The largest change in transport- 
related protein abundance was observed in anaerobically grown 
Arex cells: 8 proteins showed significant abundance level changes 
(p<0.05). This suggests that the Rex-dependent control of the 
abundance pattern of transport related-proteins was stronger in 
anaerobic fermentative conditions than in aerobic respiratory 
conditions. 

Transcriptional regulation-related proteins. Arex up-regulated 
the synthesis of 6 and 4 transcriptional regulators under 
anaerobiosis and aerobiosis, respectively. Among these, we found 
regulators of the Mer and Mar families, which control many 
ceUular processes including oxidative stress response and virulence 
[43,44]. Interestingly, we noted that ResD was significantly up- 
regulated in anaerobically grown Arex cells but remained 
unchanged in aerobicaUy grown Arex cells. ResD is involved in 
both metabolism and toxinogenesis in B. cereus [22,32]. 

Stress-response related proteins. Among the stress response- 
related proteins, the antioxidant protein, OhrA [7,9], was 
significantly up-regulated in aerobicaUy grown Arex ceUs and its 
abundance remained constant under anaerobiosis (Table 2). Like 
OhrA, the ribosome-associated heat-shock protein, Hspl5 
(BC_0062), showed, under aerobiosis, a significant change in 
abundance level. This protein could be involved in the recycling of 
free 50 S ribosomal subunits [45]. Interestingly, we noted 
significant changes in the abundance patterns of 50 S ribosomal 
subunits in aerobicaUy grown Arex ceUs (Table S6). Unlike 
aerobicaUy grown cells, anaerobically grown Arex ceUs exhibited 
higher levels of catalase, KatE (BC_1155), lower levels of the 
BC1603 protein, which is functionally annotated as a cold-shock 
protein and no significant level change of alkyl hydroperoxide 
reductase (AhpC, BC_0377). Clearly, Rex controls the abundance 
pattern of stress-related proteins in an oxygen-dependent manner. 

Phage-related proteins. Five phage-related proteins were 
strongly up-regulated in anaerobically grown Arex cells (log2> 
1.5, /'-value<0.05) and their abundance remained constant in 
aerobicaUy grown Arex ceUs. The functions of these proteins are 
unknown but they are considered to be potentiaUy beneficial to B. 
cereus in facing adverse environmental conditions [42] . 

Virulence-related proteins. Table 2 shows that rex deletion has a 
significant impact on the abundance levels of hemolysin I (Hly I), 
neutral protease Npr2 (BG_2735), protease HhoA (BC_3600), and 
SSEB (BC_1366) proteins under aerobiosis, all of which are 
predicted to be secreted virulence factors [4,46] . 

Insights into the extracellular proteome of the Arex 
mutant 

(i) Focus on toxin-related proteins. Although the intracel- 
lular abundance level of Hlyl was higher in the aerobicaUy grown 
Arex mutant strain, the extracellular level of this protein was not 

significandy different in the mutant compared with wild-type cells 
(Table S7). This was also the case for EntB, a putative enterotoxin 
[9]. Howev(;r, 12 other toxin-related proteins showed a signifi- 
candy higher abundance level in the exoproteome of aerobicaUy 
grown Arex ceUs than wUd-type ceUs (Table 3). Some of these 
proteins also showed a higher abundance level under anaerobiosis, 
but to a lesser extent, i.e. the three components of Nhe (NheA, 
NheB and NheC), EntA, EntC and the LI component of Hbl. 



From these data, we conclude that Rex regulates the toxinogenic 
profUe of the B. cereus exoproteome in response to oxygen 
availability. 

(ii) Other predicted extracellular proteins. Like the toxin- 
related proteins, many degradative enzymes, flagella components 
and proteins that are released from the cell wall, are more 
abundant in the exoproteome of the Arex mutant strain than wUd- 
type ceUs, especiaUy under aerobiosis (Table 4). AU of these 
proteins are predicted to be extraceUular proteins and, except for 
flagella components, all of them encompass signal peptides in their 
N-terminal primary structure and/or transmembrane helices. One 
metabolism-related protein, the subunit 0( of EOF 1 ATP synthase 
(AtpA), and three proteins of unknown function (BC_1894, 
BC_5360 and BC_4062) also showed a higher abundance level 
in the exoproteome of aerobicaUy grown Arex ceUs. Except for 
BC_1894, which was annotated as a prophage protein and thus 
could be secreted via holin(-like) pathway [47], these proteins 
contain a signal peptide, and thus are predicted to be secreted by 
classical export pathways. 

(iii) Predicted cytoplasmic proteins. In contrast to pro- 
teins predicted to be classically secreted, many proteins predicted 
to be cytoplasmic proteins, including metabolic proteins, stress 
response-related proteins, translation-related proteins, and other 
proteins were significantly less abundant in the exoproteome of the 
Arex mutant than the parental strain (Table 4). We noted that 
most of the cytoplasmic proteins that showed a significant 
abundance level decrease under aerobiosis did not show significant 
abundance level changes under anaerobiosis, and vice-versa 
(Table 4 and Table S7). Finally, the percentage of predicted 
cytoplasmic proteins significantly decreased in the Arex exopro- 
teome compared with the wUd-type exoproteome in the presence 
of oxygen, whUe the percentage of typically secreted proteins 
(toxin-related proteins and degradative enzymes and adhesins) was 
significandy higher (Fig. 4). 

The absence of Rex does not affect the mRNA level of 
toxin-related proteins 

An explanation for the increased abundance levels of toxin- 
related proteins in the Arex mutant exoproteome could be a 
change in their mRNA levels. We thus tested whether transcrip- 
tion of genes encoding these proteins was altered in the Arex 
mutant compared with wUd-type. Table 3 indicates that rex 
deletion significandy down-regulated entC and entFM under 
aerobiosis, but did not significandy change transcription of aU 
other toxin-related genes. These transcriptomic results are not 
concordant with the proteomic data. This suggests that toxin- 
related protein abundance may be mainly modulated at the post- 
transcriptional level when Rex is lacking. 

Discussion 

B. cereus Rex is an autlunitic Rex transcriptional factor because 
(i) it is able to bind to promoter fragments in vitro, incUcating that 
at least in some cases the effect of Rex is mediated by direct 
interaction with the promoters, (ii) it impacts the cellular 
metabolism and oxidative stress tolerance, and (iii) it acts mainly 
as a repressor at the proteomic level [12,13,14,15,20,48]. 
However, B. cereus Rex also plays an important role in fuUy 
aerobic respiratory metabolism and it modulates the toxinogenic 
profUe of exoproteome. 

When oxygen was supplied to pH-regu!ated batch cultures at 
p02 = 100%, B. cereus metabolizes glucose to carbon dioxide by 
oxidation of glycolytic pyruvate in the tricarboxylic acid (TCA) 
cycle. This reaction produces NADH, which then fuels oxidative 
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Figure 4. Functional distribution of proteins identified in Arex and wild-type B. cereus exoproteomes under aerobiosis and 
anaerobiosis. The diagrams represent the average proportion of each functional protein group based on total spectral counts (Table S5). For 
readability reasons, standard deviations (below 10%) are not shown. 
doi:1 0.1 371 /journal.pone.01 07354.g004 



phosphorylation to produce ATP, with minimal production of 
lactate and elimination of carbon excess through excretion of 
acetate. Rex deficiency facilitates the entry of carbon flow into the 
NADH-recycling lactate pathway at the expense of pyruvate 
oxidation into acetyl-CoA (Fig. 1). This promotes lactate produc- 
tion at the expense of acetate without impairment of TCA 
capacity, since acetate was always produced as a waste [28] . In the 
absence of oxygen or other external electron acceptors, ATP 
synthesis occurs at the level of substrate phosphorylation and 
NADH is reoxidized in terminal step fermentative reactions from 
pyruvate (Fig. 1). When grown in pH-controUed anaerobic batch 
cultures (p02 = 0%), B. cereus produces lactate as the main 
glucose by-product to satisfy the demand for redox balance. Rex 
deficiency favored the production of more reduced metabolites, as 
evidenced by the increase in the ethanol to acetate ratio and the 
decrease of lactate production. This indicates that Rex regulates 
the carbon flow distribution at the pyruvate node by favoring 
carbon flow through the NADH-recycling lactate pathway at the 
expense of Pfl-dependent fermentative pathways under anaerobi- 
osis. A key question is how Rex controls the carbon flow into the 
NADH-recycling lactate pathway in B. cereus cells grown in fuUy 
aerobic respiratory and anaerobic fermentative conditions? This 
control process could involve interplay of transcriptional and/ or 
post-transcriptional regulation of key glycolytic enzymes (such as 
Pgm, BC_4919 under anaerobiosis and PK, BC_4599 under 



aerobiosis) and post-translational regulation of lactate dehydroge- 
nase. Evidendy, further intensive work is required to unravel the 
complexity of the regulatory network involving Rex. Finally, by 
controlling the entry of carbon flow into the lactate pathway. Rex 
controls the availability of glycolytic intermediates for macromo- 
lecular synthesis as well as supporting NADPH production under 
both aerobiosis and anaerobiosis. 

In contrast to anaerobic fermentation, aerobic respiration is a 
major source of reactive oxygen species (ROS) generation [49]. B. 
cereus, like other facultative anaerobes, uses scavenging systems to 
control the level of ROS [50]. Sequestration of ROS leads to 
oxidative deactivation of these scavenging systems, which are then 
reactivated direcdy or indirectly depending on NADPH level [51]. 
By modulating the distribution pattern of antioxidant enzymes 
(Tables 2 and 4), and possibly NADPH production. Rex may thus 
modulate the effectiveness of antioxidative defense systems. 
Specifically, Rex may maximize antioxidant system activity under 
anaerobiosis while restricting this activity under aerobiosis; this 
could explain why the resistance of fi. cereus cells towards external 
H2O2, when Rex is absent, is lower under anaerobiosis and higher 
under aerobiosis. 

Except for Hlyl, the abundance levels of classical extracellular 
proteins (such as toxins) were enhanced in B. cereus exoproteome 
when Rex was absent, especially under aerobiosis. In B. cereus, as 
in other Gram-positive bacteria, most extracellular proteins are 
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synthesized as precursors, with typical N-terminal signal peptides, 
and exported from the cytoplasm by the Sec-dependent pathway 
[52]. This is the case for toxins [6]. However, unlike other toxins, 
Hlyl contains a putative N-terminal twin arginine motif that could 
target the protein to the Tat secretion pathway instead of the Sec- 
dependent pathway [53,54]. Besides classical export routes, B. 
cereus possesses speciaUzed protein export systems such the 
flagellar type III system [5,55], which translocates most flagellar 
components across the cytoplasmic membrane, and possibly the 
ESAT-6 secretion system, which mediates the secretion of virulence 
factors belonging to the WXGIOO protein family, such as the 
protein BC_5360 [56,57]. A feature of these two specialized 
secretion systems is the lack of a classical signal peptide to direct the 
substrate protein for secretion. However, like the Sec-dependent 
pathway and unlike the TAT-dependent pathway, these protein 
export systems require chaperones to prevent premature folding, 
aggregation and degradation by cytoplasmic proteases [58,59]. Rex 
may thus influence protein-assisted export by modulating the 
distribution pattern of the components of export systems as 
observed in our proteomic data (Table 2, [60]). In addition, it has 
been reported that chaperones may be targets for intracellular ROS 
[61]. Keeping the functional integrity of these chaperones, 
specifically in the context of secretion, may thus depend on Rex- 
dependent NADPH availability. The role and contribution of 
NADPH in classical and specialized secretory mechanisms remains 
an open question that deserves further experiments. 

Our results show that many cytoplasmic proteins predicted to 
be imoh'cd in active growth, either directly (such as metabolic 
enzymes, translational and post-translational-related proteins) or 
indirectly by preventing inappropriate gene expression (such as 
AbrB and CodY under aerobiosis [62]), were less abundant in B. 
cereus exoproteome when Rex was absent. We assume that cell 
lysis cannot account for changes in the abundance level of these 
cytoplasmic proteins since: (i) cell viability count did not change 
(data not shown), (ii) EF-Tu did not accumulate in extracellular 
medium (Table S7, [63]), and (iii) the abundance pattern of 
cytoplasmic proteins was dependent on growth conditions and did 
not always correlate with intracellular abundance changes 
(Table 2 and 4). Among the cytoplasmic proteins identified in B. 
cereus exoproteome, several have already been identified in the 
exoproteome of other bacteria, and many have been described as 
extracellular moonlighting components playing a role in bacterial 
virulence [64,65] . The non classical mechanism of secretion that 
could explain their release remains to be firmly described [47,66]. 
One hypothesis is that signal-less predicted cytoplasmic proteins 
might be secreted within outer membrane vesicles [63]. Such 
nanovesicles have been observed in Bacillus anthracis [67]. 
Interestingly, ATP is the primary metabolic factor involved in 
secretory vesicle movements. ATP generation is lower in 
fermenting cells than in non-fermenting cells, according to 
metabolism. In addition. Rex promotes ATP generation by 
enhancing NADH-producing pathways at the expense of NADPH 
pathways. Therefore, if ATP availability is a triggering factor for 
export of cytoplasmic proteins, this could explain the higher export 
of cytoplasmic proteins in aerobicaUy grown cells compared to 
anaerobicaUy grown cells, as well as the substantial decrease of 
cytoplasmic protein export under aerobiosis when Rex was absent. 

In conclusion, Rex in B. cereus plays a pivotal role in controlling 
the cross-talk between the metabolic networks that control growth, 
oxidative defense machinery and extracellular accumulation of 
toxins. In the context of B. cereus natural environment. Rex may 
thus help the pathogen to (i) maintain its host alive until 
transmission to the next host can be achieved, by controlling its 
growth and production of its toxins [68], and (ii) better survive 



stress conditions under anaerobiosis, albeit to the detriment of 
maximizing its survival when oxygen is present [10,69,70]. 

Supporting Information 

Figure SI Genetic organization of rex chromosomal 
region of B. cereus ATCC 14579. Panel A. Large arrows 
represent open reading frames identified in strain ATCC 14579, 
and their orientation shows the transcriptional direction. Panel B. 
Nucleotide sequence of rex ORF (black capital letters), as well as 
promoter and terminator regions (blue letters), are shown. 
Transcriptional start site (+1), putative —35 and —10 motifs and 
putative terminator are underlined. Start and stop translation 
codons are in bold. 
(TIF) 

Figure S2 Survival of B. cereus F4430/73 cells towards 
external hydrogen peroxide. Cells were grown in liquid 
cultures to mid-exponential growth phase either aerobicaUy or 
anaerobicaUy, and subjected to 20 and 10 mM HjOj stress, 
respectively. Samples were taken at 0 and 20 min. Colony forming 
units per mL were counted and expressed as (N/No)xl00. Error 
bars represented the standard deviation from three independent 
measures. Significant differences (p<0.05) between mutant and 
wild-type strains are indicated with asterisks. 
(TIF) 

Figure S3 SDS-PAGE analysis of overproduced and 
purified B. cereus Rex. Rex purification fractions were 
analyzed by electrophoresis on an SDS-12% polyacrylamide gel 
after Coomassie Brillant Blue staining (Panel A) or Imperial 
Protein Stain (Thermo scientific, Panel B). Position and molecular 
weights (kDa) of markers (lanes 1) are given on the left. Panel A: 
Lane 1, standard proteins; Lane 2, soluble whole cell extract from 
E. coli. Panel B: Lane 1, standard proteins; Lane 2, purified 
protein after Co^"^ IMAC. 
(TIF) 

Figure S4 Venn diagrams summarizing proteins that 
showed significant abundance level changes in the 
cellular and extracellular proteomes of \rex compared 
with wild-type B. cereus cells under anaerobiosis (blue) 
and aerobiosis (red). Proteins up- and down-regulated by Arex 
were distinguished. 
(TIF) 

Table SI Results from controlled batch cultures ot Arex 
mutants and its parent strain, B. cereus F4430/73. 

(DOCX) 

Table S2 List of assigned MS/MS spectra from the 
cellular proteomes. 

(XLSX) 

Table S3 List of assigned MS/MS spectra from the 
exoproteomes . 

PCLSX) 

Table S4 List of proteins identified in cellular pro- 
teomes and their normalized spectral count. 

PCLSX) 

Table S5 List of proteins identified in exoproteomes 
and their normalized spectral count. 

PCLSX) 

Table S6 Fold changes of cellular proteins and their 
statistical significance. 

pCLSX) 
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Table $7 Fold changes of exoproteins and their statis- 
tical significance. 

(XLSX) 

Acknowledgments 

S.L. held a fellowship from the Institut National de la Recherche 
Agronomique and the French Region of Provence Alpes Cotes d'Azur 

References 

1. Bottonc EJ (2010) Bacillus cereus, a volatile human pathogen. Clin Microbiol 

Rev 23: 382-398. 

2. Stcnfors Arncsrn LP. Fagcrlund A. Ciranum PE (2008) Erom soil to gut: Bacillus 
cereus and its food poisoning toxins. EEJMS Microbiol Rev 32: 579—606. 

3. Ccuppcns S, Boon N, Uyttrndacic JM (2013) Divrrsity of Bacillus corns group 
strains is rellcctcd in thrir broad range of pathogenicity and diverse ecological 
lifestyles. FEMS Microbiol Ecol. 84: 433-50. 

4. Ramarao N, Sanchis V (20 1 3) The pore-forming haemolysins of Bacillus cereus: 
a review. Toxins (Basel) 5: 1119—1139. 

5. Senesi S, Ghelardi E (2010) Production, secretion and biological activity of 
Bacillus cereus entcrotoxins. Toxins (Basel) 2: 1690-1703. 

6. Fagerlund A, Lindback 'E, (iranum PE (2010) Bacillus cereus cytotoxins Hbl, 
Nhc and CvtK are secreted \'ia the Sec translocation pathwav. BMC Microbiol 
10: 304. 

7. Clair G, Eorphelin A, Armengaud J, Duport C (2013) OhrRA functions as a 
rcdox-responsi\e system controlling toxinogenesis in Bacillus cereus. J Proteo- 
mics 94: 527-539. 

8. Clair G, Roussi S, Armengaud J, Duport C (2010) Expanding the known 
repertoire of virulence factors produced by Bacillus cereus through early 
secretome profiling in three redox conditions. Mol Cell Proteomics 9: 1486— 
1498. 

9. Clair G, Armengaud J, Duport C (2012) Restricting Fermentative Potential by 
Proteome Remodeling: An adaptative strategy evidenced \n Bacillus cereus. Mol 
Cell Proteomics 11: Mill 013102. 

10. Marteyn B, West NP, Browning DE, ColeJA, ShawJG, et al. (2010) Modulation 
of Shigella virulence in response to available oxvgen in \'ivo. Nature 465: 355— 
358. 

11. Moriarty-Craige SE, Jones DP (2004) Extracellular thiols and thiol/disuUide 
redox in metabolism. Annu Rev Nutr 24: 481-509. 

12. Larsson J'E, Rogstam A, von Wachenfeldt C (2005) Coordinated patterns of 
cytochrome bd and lactate dehydrogenase expression in Bacillus subtilis. 
Microbiology 151: 3323-3335. 

13. Brekasis D, Paget MS (2003) A novel sensor of NADH/NAD+ redox poise in 
Streptomyces coelicolor A3(2). EMBO J 22: 4856-4865. 

14. Gyan S, Shiohira Y, Sato I, Takeuehi M, Sato T (2006) Regulatopy- loop 
between redox sensing of the NADH/1VAD(+) ratio bv Rex (YdiH) and 
oxidation of NADH by NADH dehydrogenase Ndh in Bacillus subtilis. 
J Barteriol 188: 7062-7071. 

15. Pagcls M, Euchs S, Pane-Farre J, Kohler C, Mcnschner L, et al. (2010) Redox 
sensing by a Rex-family repressor is involved in the regulation of anaerobic gene 
expression in Staphylococcus aureus. Mol Microbiol 76: 1142-1161. 

16. Sickmier EA, Brekasis D, Paranawithana S, BonannoJB, Paget MS, et al. (2005) 
X-ray structure of a Rex-family repressor/NADH complex insights into the 
mechanism of redox sensing. Structure 13: 43—54. 

17. Somerville GA, Proctor RA (2009) At the crossroads of bacterial metabolism and 
virulence factor synthesis in Staphylococci. Microbiol Mol Biol Rev 73: 233-248. 

18. V'esic D, Kristich CJ (2013) A Rex familv transcriptional repressor influences 
H202 accumulation hy Enterococcus faeca lis. j Bacteriol 195: 1815-1824. 

19. Wang E, Ikonen TP, Knaapila M, Svergun D, Logan DT, et al. (2011) Small- 
angle X-ray scattering study of a Rex family repressor: conformational response 
to NADH and NAD+ binding in solution. J Mol Biol 408: 670-683. 

20. Bitoun JP, Nguyen AH, Fan Y, Burne RA, Wen ZT (2011) Transcriptional 
repressor Rex is involved in regulation of oxidative stress response and biofilm 
formation by Streptococcus mutans. FEMS Microbiol Lett 320: 110-117. 

21. Ravcheev DA, Li X, Latif H, Zengler K, Leyn SA, et al. (2012) Transcriptional 
regulation of central carbon and energy metabolism in bacteria by redox- 
responsivc repressor Rex. J Bacteriol 194: 1145—1157. 

22. EsbelinJ, ArmengaudJ, Zigha A. Duport C (2009) ResDE-dependent regulation 
of enterotoxin gene expression in Bacillus cereus: evidence for multiple modes of 
binding for ResD and interaction with Fnr. J Bacteriol 191: 4419-4426. 

23. Esbelin J, Jouanncau Y, ArmengaudJ, Duport C (2008) ApoEnr binds as a 
monomer to promoters regulating the expression of enterotoxin genes of Bacillus 
c^reM5. J Bacteriol 190: 4242^251. 

24. Zigha A, Rosenfeld E, Schmitt P, Duport C (2007) The redox regulator Fnr is 
required for fermentative growth and enterotoxin synthesis in Bacillus cereus 
F4430/73. J Bacteriol 189: 2813-2824. 

25. Laouami S, Mcssaoudi K, Alberto E, Clavel T, Duport C (2011) Lactate 
dehydrogenase A promotes communication between carbohydrate catabolism 
and virulence in Bacillus cereus.} Bacteriol 193: 1757-1766. 



and G.C. held a fellowship from the Ministere de la Recherche et de 
TEnseignement Superieur. 

Author Contributions 

Conceived and designed the experiments: SL JA CD. Performed the 
experiments: SL GC CD. Analyzed the data: SL JA CD. Contributed 
reagents/materials/analysis tools; SL GC JA CD. Contributed to the 
writing of the manuscript: CD. 



26. Esbehn J, Jouanncau Y, Duport C (2012) Bacillus cereus Enr binds a [4Fe-4 S] 
cluster and forms a ternary complex with ResD and PlcR. BMC Microbiol 12: 
125. 

27. Ivanova N, Sorokin A, Anderson I, Galleron N, Candelon B, et al. (2003) 
Genome sequence of Bacillus cereus and comparative analysis with Bacillus 
anthracis. Nature 423: 87-91. 

28. Rosenfeld E, Duport C, Zigha A, Schmitt P (2005) Characterization of aerobic 
and anaerobic vegetative growth of the food-borne pathogen Bacillus cereus 
F4430/73 strain. CanJ Microbiol 51: 149-158. 

29. Murphy E (1985) Nucleotide scf^uence of a spectinomycin adenyltransferase 
AAD(9) determinant from Staphylococcus aureus and its relationship to AAD(3") 
(9). Mol Gen Cenet 200: 33-39. 

30. Arnaud M, Chastanet A, Dcbarbouille M (2004) New vector for elficient allelic 
replacement in naturally nontransformable, low-GC -content, gram-positive 
bacteria. Appl Environ Microbiol 70: 6887-6891. 

31. Arantes O, Lereclus D (1991) Construction of cloning vectors for Bacillus 
thuringiensis . Gene 108: 115-119. 

32. Duport C, Zigha A, Rosenfeld E. Schmitt P (2006) Control of enterotoxin gene 
expression in Bacillus cereus El 1311/73 involves the redox-sensitive ResDE 
signal transduction system. J Bacteriol 188: 6640 6651. 

33. Car\/alho PC:, Eischer jS^ Chen El, Yairs JR, 3rd. Barbosa VC (2008) 
PatternLab for proteomics: a tool for dili'ereniial shotgun proteomics. BMC 
Bioinformatics 9: 316. 

34. Vizcaino JA, Cote RG, Gsordas A, Dianes JA, Eabregat A, et al. (2013) The 
PRoteomics IDEntifications (PRLDE) database and associated tools: status in 
2013. Nucleic Acids Res 41: D1063-1069. 

35. Wang Q, Ou MS, Kim Y, Ingram LO, Shanmugam KT (2010) Metabolic flux 
control at the pyruvate node in an anaerobic Escherichia coli strain with an 
active pyruvate dehydrogenase. Appl Environ Microbiol 76: 2107-2114. 

36. Wang E, Bauer MC, Rogstam A, Linse S, Logan DT, et al. (2008) Structure and 
functional properties of the Bacillus subtilis transcriptional repressor Rex. Mol 
Microbiol 69: 466-478. 

37. Schau M, Chen Y, Hulett EM (2004) Bacillus subtilis YdiH is a direct negative 
regulator of the cyd^\BCD operon. J Bacteriol 186: 4585-4595. 

38. Sprynski N, Eelix C, O'Callaghan D, Vergunst AC (20 1 2) Restoring virulence to 
mutants lacking subunits of multiprotein machines: functional complementation 
of a. Brucella virB5 mutant. FEBS Open Bio 2: 71-75. 

39. Spira WM, Goepfert JM (1975) Biological characteristics of an enterotoxin 
produced hy Bacillus cereus. CanJ Microbiol 21: 1236- 1246. 

40. Garrigues C, Goupil-Eeuillerat N, Cocaign-Bousc^uct M, Renault P, Lindley 
ND, et al. (2001) (ilucose metabolism and regulation of glvcolysis in Lactococcus 
lactis strains with decreased lactate dehydrogenase activity. Metab Eng 3: 211— 
217. 

41. Wichgers Schreur PJ, Rebel JM, Smits MA, van PuttenJP, Smitii HE (201 1) Lgt 
processing is an essential step in Streptococcus suis lipoprotein mediated innate 
immune activation. PLoS One 6: e22299. 

42. Wang X, Kim Y, Ma Q, Hong SH, Pokusaeva K, et al. (2010) Cryptic 
prophages help bacteria cope with adverse environments. Nat Gommun 1: 147. 

43. HoUey TA, Stevenson CE, Bibb MJ, Lawson DM (2013) High resolution crystal 
structure of Sco5413, a widespread actinomycete MarR family transcriptional 
regulator of unknown function. Proteins 81: 176—182. 

44. Newberrv KJ, Brennan R(i (2004) Ihe structural mechanism for transcription 
activation by McrR family member multidrug transporter activation, N 
terminus. J Biol Chem 279: 20356-20362. 

45. Jiang L, Schaffitzel C, Bingel-Erlenmeyer R, Ban N, Korber P, et al. (2009) 
Recycling of aborted ribosomal 50 S subunit-nascent chain-tRNA complexes by 
the heat shock protein Hspl5. J Mol Biol 386: 1357-1367. 

46. Gohar M, Facgri K, Perchat S, Ravnum S, Okstad OA, et al. (2008) The PlcR 
virulence regulon of Bacillus cereus. PLoS One 3: e2793. 

47. Tjalsma H, Antelmann H, Jongbloed JD, Braun PG, Darmon E, et al. (2004) 
Proteomics of protein secretion hy Bacillus subtilis: separating the "secrets" of 
the secretome. Microbiol Mol Biol Rev 68: 207-233. 

48. Bitoun JP, Liao S, Yao X, Xie GG, Wen ZT (2012) The redox-sensing regulator 
Rex modulates central carbon metaboHsm, stress tolerance response and biofilm 
formation by Streptococcus mutans. PLoS One 7: e44766. 

49. Seaver LC, Imlay JA (2001) Alkyl hydroperoxide reductase is the primary 
scavenger of endogenous hydrogen peroxide in Escherichia coli. J Bacteriol 183: 
7173-7181. 

50. Mols M, Abee T (2011) Primary and secondary oxidative stress in Bacillus. 
Environ Microbiol. 13: 1387-1394. 



PLOS ONE I www.plosone.org 



16 



September 2014 | Volume 9 | Issue 9 | e107354 



Proteomics of B. cereus Arex 



51. MailloLix RJ, Lcmirc J, Appanna \ D (2011) Metabolic networks to combat 
oxidative stress in Pseudomonas jluorescens. Antonie Van Lccuwenhoek 99: 
433-442. 

52. Schneewind O, Missiakas D (2013) Sec-secretion and sortase-mediated 
£inchoring of proteins in Gram-positive bacteria. Biochim Biophys Acta. 

53. Yuan J, Zweers JC, van DijlJM, Dalbey RE (2010) Protein transport across and 
into cell membranes in bacteria and archaea. Cell Mol Life Sci 67: 179-199. 

54. Goosens VJ, Monteferrante CG, van DijlJM (2013) The Tat system of Gram- 
positive bacteria. Biochim Biophys Acta. 

55. Senesi S, Salvetti S, Celandroni F, Ghelardi E (2010) Features of Bacillus cereus 
swarm cells. Res Microbiol 161: 74!?-749. 

56. Burts ML, Williams VVA, DeBord K, Missiakas DM (2005) EsxA and EsxB arc 
secreted by an ESAT-6-like system that is required for the pathogenesis of 
Staphylococcus aureus infections. Proc Natl .Vcad Sci U S A 102: 1169—1174. 

57. Fallen MJ (2002) The ESAT-6/\VXG100 superfamily - and a new Gram- 
positive secretion system? Trends Microbiol 10: 209—212. 

58. Alonzo F, 3rd, Xayarath B, WhisstockJC, Frcitag NE (201 1) Functional analysis 
of the Listeria monocytogenes secretion chaperone PrsA2 and its multiple 
contributions to bacterial virulence. Mol Microbiol 80: 1530-1548. 

59. Moran-Barrio J, Limansky AS, Viale AM (2009) Secretion of GOB metallo-beta- 
lactamase in Escherichia coli depends strictly on the cooperation between the 
cytoplasmic DnaK chaperone system and the Sec machinery: completion of 
folding and Zn(ll) ion acquisition occur in the bacterial periplasm. Antimicrob 
Agents Chemother 53; 2908-2917. 

60. Roberts IN, Lam XT, Miranda H, Kieselbach T, Funk C (2012) Degradation of 
PsbO by the Deg protease HhoA Is thioredoxin dependent. PLoS One 7: 
e45713. 



61. Ezraty B, Aussel L, Barras F (2005) Methionine sulfoxide reductases in 
prokaryotes. Biochim Biophys Acta 1703: 221-229. 

62. Frenzel E, Doll V, Pauthner M, Lucking G, Scherer S, et al. (2012) CodY 
orchestrates the expression of virulence determinants in emetic Bacillus cereus 
by impacting key regulatory circuits. Mol Microbiol 85: 67—88. 

63. Meneses N, Alendoza-Hernandez G, Encamacion S (2010) The extracellular 
proteomc oi' lihizohi mil etli CE3 in exponential and stationary growth phase. 
Proteome Sci 8: 51. 

64. Copley SD (2012) Moonlighling is mainstream; paradigm adjuslmcnl required. 
Biocssays 34: 578-588. 

65. Henderson B, Martin A (2013) Bacterial moonlighting proteins and bacterial 
virulence. Curr Top Microbiol Immunol 358; 155-213. 

66. Yang CK, Ewis HE, Zhang X, Lu CD, Hu HJ, et al. (2011) Nonclassical protein 
secretion by Bacillus siihtilis in the stationary phase is not due to cell lysis. 
J Bacteriol 193; 5607-5615. 

67. Rivera J, Cordero RJ, Nakouzi AS, Frases S, Nicola A, et al. (2010) Bacillus 
anthracis produces membrane -derived vesicles containing biologically active 
toxins. Proc Natl Acad Sci USA 107; 19002-19007. 

68. ten Bokum AM, Movahedzadeh F, Frita R, Bancroft GJ, Stoker NG (2008) The 
case for hypervirulence through gene deletion in Mycobacterium tuberculosis. 
Trends Microbiol 16: 436-441. 

69. Jakesevic M, Aaby K, Borge GI, Jeppsson B, Ahme S, et al. (201 1) Antioxidative 
protection of dietary bilberry, chokeberry and Lactobacillus plantarum HEAL19 
in mice subjected to intestinal oxidative stress by ischemia-reperfusion. BMC 
Complement i\ltern Med 11:8. 

70. Marteyn B, Seorza FB, Sansonetti PJ, 'Fang C (2010) Breathing life into 
pathogens: the influence of oxygen on bacterial virulence and host responses in 
the gastrointestinal tract. Cell Microbiol 13: 171-176. 



PLOS ONE I www.plosone.org 



17 



September 2014 | Volume 9 | Issue 9 | e107354 



